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Abstract 

We investigated the dependence of the K~d scattering length on models of KN interaction with 
one or two poles for A(1405) resonance. The KNN — ttTjN system is described by coupled-channel 
Faddeev equations in AGS form. Our new two-body KN — 7rS potentials reproduce all existing 
experimental data on K~p scattering and kaonic hydrogen atom characteristics. New models of 
T,N — AA^ interaction were also constructed. Comparison with several approximations, usually 
used for scattering length calculations, was performed. 

PACS numbers: 13.75.Jz, ll.80.Gw 
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I. INTRODUCTION 



Investigation of K~d system can shed more light on the KN interaction, necessary for 
study of antikaonic nuclear clusters, which attracted large interest recently [l| • The interac- 
tion is not very well known, in particular, there are debates about the nature of the A (1405) 
resonance. The question is whether it is a single resonance in 7rE and a quasi-bound state 
in KN channel or the bump, which is usually understood as A(1405) resonance, is an effect 
of two poles. The advantage of the K~d system is the possibility of proper description of 
its dynamics using Faddeev equations j^J. 

Recently we constructed coupled-channel KN — 7rS potentials in one- and two-pole 
form [3|, which reproduce all existing experimental data on K~p scattering and K~p atom 
characteristics equally well, therefore it is not possible to give preference to any of the ver- 
sions. A possible way to clarify the question concerning the nature of the A(1405) resonance 
is to perform few- or many-body calculations using one- and two-pole KN — 7rS poten- 
tials as an input. Having in mind SIDDHARTA experiment [4j, measuring characteristics 
of kaonic deuterium atom, we calculated K~d scattering length a K - d and investigated the 
dependence of the results on the models of KN interaction with newly obtained parameters. 
The scattering length gives possibility to calculate kaonic deuterium level shift and width. 
Comparison of the theoretical results with experimental ones could allow to choose between 
the two KN — 7rS interaction versions. 

Dependence of ax-d on other two-body interactions, necessary for the description of 
the KNN — ttTjN system, was also investigated: we used several models of A^A" (with 
and without short range repulsion) and HN(— AN) interactions. In addition to the full 
coupled-channel calculation we performed checks of commonly used approximations for K~d 
scattering. In particular, we solved one-channel Faddeev equations using exact optical and 
simple complex KN potentials approximating the KN — 7rS models of interaction. We also 
checked the "Fixed center approximation to Faddeev equations" formula. 

The formalism used for the coupled-channel K~d scattering length calculation is described 
in the next section. Section HTT1 is devoted to the two-body input: the description of the one- 
and two-pole KN — 7rS potentials with newly obtained parameters in the first subsection is 
supplemented with additional arguments for equivalence of the two versions. The following 
subsections of Section llHl are devoted to A^A" and SA^ — AA" potentials. Section [TVl contains 
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information about approximate methods, usually used in K~d scattering length calculations. 
The full and approximate results are shown and discussed in Section [V] while Section [VI] 
concludes the paper. 



II. COUPLED-CHANNEL AGS EQUATIONS FOR KNN - vrSiV SYSTEM 

As in 6|,|7| we directly include nUN channel into original three-body Faddeev equations in 
the Alt-Grassberger-Sandhas (AGS) form [5], which leads to the coupled-channel equations: 

U?f = <W (i - k) (Go)" 1 + E (1 - Su.) IT Gl Ug, (1) 

fe,7=l 

where "particle channel" indices a, (3 = 1,2,3 are introduced in addition to the usual Faddeev 
partition indices i,j, k = 1, 2, 3, see Table I of [7|. The equations define unknown operators 
lHj* , describing the elastic and re-arrangement processes j 13 + (k^i 13 ) — > i a + (j a k a ). The free 
Green's function is diagonal in channel indices: G^ 13 = 5 a/ 3 Gq. The inputs for the system of 
equations {TJ are two-body T-matrices, embedded into three-body space: describes the 
interaction between the particles j and k (i ^ j ^ k) in channels a, (3. Like in jjj], here we 
have T t NN , T[ N and Tp N , which are usual one-channel two-body T-matrices in three-body 
space, describing NN, irN, and EiV interactions, respectively. The T t KK , T/ r7r , Tf K , and 
T, K7T are elements of the coupled-channel T-matrix for the KN — 7rE system. 

Q 

In contrast to the calculation of the quasi-bound K pp state [7[ , where one-term isospin 
(/) dependent separable potentials were used, now we write AGS equations for N-term 
isospin dependent separable potentials 

ATf 

V iJ = ^2 K(m),I \9i(m),l)(9i(m),l\> ( 2 ) 
m=l 

which lead to separable T-matrices 

ATf 

^if = ^2 \9i{m),l) r i^nn),l{di(n),l\ ■ (3) 
m,n=l 

Here N" is a number of terms of the separable potential, A is a strength constant, while g 
is a form-factor. Bound state wave function of the two-body subsystem, described by such 
a potential, has the form 

Wi,l) = C i{m),I G o( z bnd) \9i{m),l) > ( 4 ) 

m=l 



where the coefficients C9is j are constants and Zb n d is a binding energy. We used two slightly 
different versions of a two-term nucleon-nucleon separable potential (iV" = 2) in the ax-d 
calculations. All other models of interactions, one more V NN among them, are one-term 
potentials with N° = 1. 

The amplitude f^j.j. of the K(NN) — » K(NN) reaction with NN in isospin-zero state 
in the initial and final states is defined by the following matrix element: 

fu,00(Pl,Phztot) = -(27T)Vi (ph^lolUllooiZtatM^p 1 !) , (5) 

where [if is the three-body reduced mass defined by 

mf{m a + ml) . . 

mf + m" + m£ 

Pi and p' x are initial and final relative momenta of antikaon with respect to the NN pair, 
correspondingly, while z tot is the total energy of the three-body system. Deuteron wave 
function ip\ in Eq.(j5]) is defined by Eq.(jl]), the transition operator Uff^j. has two additional 
isospin indices compared to Eq.fJTJ. The scattering length is the amplitude at zero kinetic 
energy z% in : 

a K -d = fufloiPi -> 0,Pi -> 0; z tat -» z] h ) (7) 

with z° h = X]^=i m ? being the KNN (a = 1) or nY,N (a = 2) threshold energy (the total 
energy is defined by z tot = z% in + zf h ). 

Taking into account forms of T-matrices ([3]) and wave functions (j4j), introducing new 
operators 



rap _ i a \na jjap i , B 

i{m)j,hlj - WiM./J'^O ' /.;./ /, ' /./ 

Z i(m)j{n),hl 3 = <W ^j(m)j(n),; ( / j = <W i 1 ~ (9i{m),h\ G \9*j{n),I 3 ) ) ( 9 ) 
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and substituting them into the system flTJ, we can write a system of equations for the new 
unknown operators X^ LI ,\ 



Nf 



m=l 

3 

+ ^i(l)k(m),lil k T k(mn),I k X k(n)j,I k Ij ' 



fe,7=l m,n=l / fe 
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The number of equations in the system ( FlUj) is defined by the number of all form-factors g. 
Therefore, the system fflUj) with a two-term NN and one-term other potentials consists of 
20 equations. 

Two identical nucleons, entering the first (KNN) channel, require antisymmetrization of 
the system of equations. Orbital momentum of all two-body interactions was set to zero. 
The main KN — 7r£ potential was constructed with orbital momentum / = since the 
interaction is dominated by the s-wave A(1405) resonance. The interaction of 7r-meson with 



the nucleon is mainly in p-wave, however, as it was shown in [38], the addition of "small two- 
body interactions", nN among them, changes the resulting a^-d very slightly (of the order 
of 1% or even less, see Table XIII of the paper). On the other hand, s-wave nN interaction 
is even weaker, therefore, we omitted nN interaction in our equations. Information about 
UN interaction is very poor, and there is no reason to assume significant effect of higher 
partial waves. Finally, NN interaction was also taken in Z = state only since we do not 
see physical reasons for sufficient effect of higher partial waves in the present calculation. 

Antisymmetric s-wave deuteron wave function has zero isospin and spin equal to one. 
Due to this K~d system, in contrast to K~pp , has total three-body spin (and total 
momentum) equal to one, while both KNN systems have total isospin I = \- Therefore, in 
the K~d case antisymmetrization leads to the following new operators: 

yl,asm y\ yl,asm yi __ y\ 

A l(m),0 ~~ A l(m),0' A 2,l ~~ A 2,J 3,1; 

v 2,asm _ Y 2 i v'i y2,asm _ v 2 y3 

A 1 I " A ii+ A i I) A i 3 — A i3 — A i3, 

L '2 >2 >2 1 '2 >2 >2 

-y-2,asm y2 y?j -y2,asm y2 yS 

A 2,/ —^2,1 A 3,/' A 3,/ — A 3,7 A 2,i" • 

It is necessary to note, that the K°nn state drops out from the system of equations ffTQ]) 
after the antisymmetrization because the two neutrons are in isospin one state. Therefore 
the K°nn channel has another value of the three-body spin (S = 0) than K~d (S = 1) or 
the neutrons does not satisfy Pauli principle. 

Finally, K~d scattering length can be found from 

2 

a K - d = -(2tt) 2 fjL\ Cl (m)fi X;$J(0, 0; z] h ) . (12) 

m=l 

The operator system ffTOl) written in momentum space turns into a system of integral equa- 
tions. In order to solve the inhomogeneous system we transformed the integral equations 
into algebraic ones. It is known (see. e.g. js]), that integral Faddeev equations have mov- 
ing logarithmic singularities in the kernels when scattering above a three-body threshold 



(zkin > 0) is described. "Usual" (one-channel) scattering length calculations are free of 
the singularities. For the KNN — 7rEiV system, however, where z^ h < z]: h , the permanently 
opened 7rEiV channel causes appearance of logarithmic singularities in K~d scattering length 
calculations. In the numerical procedure we handle them using the method suggested in [9(. 
The main idea of the method consists in interpolating the unknown solutions (in the interval 
containing the singular points) by certain polynomials and subsequent analytic integration 
of the singular part of the kernels. 



III. TWO-BODY INPUT 

The separable KN — 7rS, NN, and XA potentials (T5]) in momentum representation have 
the form: 

V**{k*, k>*) = £ Ag y gf m)}I (k*) gf m)>1 (^). (13) 

m=l 

Here for convenience new indices of two-body channels a,/3 = 1,2 were introduced. Corre- 
spondence between a two-body index a and a pair (a, i) of three-body channel and Faddeev 
indices, defining an interacting pair, can be established with the help of Table I of [7]. As 
before, N a defines the number of terms of the potential. As was already stated, we neglected 
here the nN interaction due to its smallness. 

In addition to the coupled-channel potentials (Eq. ffl3|) with a, j3 > 1) we used optical 
and complex one-channel potentials corresponding to them. Having in mind, that nowadays 
many authors misuse the term "optical" to a complex potential, we will call our one-channel 
potentials "exact optical" and "simple complex". An exact optical potential by defini- 
tion reproduces the elastic part of the coupled-channel interaction exactly. In particular, 
the imaginary part of the corresponding amplitude becomes zero below the lowest channel 
threshold. 

Exact optical one-channel potential, corresponding to a two-channel V (with iV a = 1), is 
given by equation (TlBl) with a, f3 — 1 and the strength parameter defined as 

A ii,o P t = A n (A) 2 ) 2 (^ 2 V 2) M) (u) 
1 l-Af<^|G«(z«)l^>' 

where Aj'^ are strength parameters of the two-channel potential, \gj) is the form- factor of 

(2) 

the second channel. Having in mind, that two-body free Green's function Gq depends on 
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the corresponding two-body energy z^ 2 \ we see, that A} 1,0pt of the exact optical potential is 
an energy-dependent complex function. In contrast to it, a strength parameter A} 1 ' omplex of 
a simple complex potential is a complex constant, therefore, the simple complex potential 
is energy independent. The strength parameter of a simple complex potential is chosen in 
such a way, that the potential reproduces some characteristics of the full interaction, say, 
scattering lengths. Exact optical and simple complex potentials take into account flux losses 
into inelastic channels through imaginary parts of the strength parameters. 



A. KN - 7r£ potential 

Two versions of phenomenological coupled-channel KN—irY, potential, constructed in J^j], 
have one- and two-pole form of A(1405) resonance and simultaneously reproduce all existing 
experimental data. In the present work we performed new fits of the experimental data 
with the same potential forms. The one-term (iV a = 1) two-channel potential is defined 
by Eq. ( TT3|) . where a = 1 denotes KN, a = 2 - ttT, channel. All physical values for 
data fitting were obtained by solving coupled-channel Lippmann-Schwinger equations with 
direct inclusion of the Coulomb potential into the K~p system. Another source of isospin 
symmetry breaking is the use of physical masses for K~, K°,p and n. 

In comparison to [3], where one one-pole and one two-pole potentials were constructed, 
here we obtained two sets of potential parameters A"^ and fif: one set for one-pole and 
another set for two-pole structure of A (1405) resonance. Each potential of these sets gives 



medium value for threshold branching ratios 10j, III] : 



7 = 2.36 ±0.04, (15) 
R n x = 0.709 ±0.011, (16) 

where the second is a new ratio, constructed from experimentally measured R c and R n : 

^ _ ^£ fyj^ 

1 — R n (1 — R c ) 

In contrast to R c and R n , the new branching ratio R^ does not contain the 7r°A channel, 
which is taken into account in our formalism only effectively through the non-zero imaginary 
part of \f K parameter. 

Elastic and inelastic K~p cross-sections K~p — > K~p, K^p — > K°n, K~p — > 7r + £~, 
K~p — ► 7r~£ + , and K~p — >■ 7r°S° are properly reproduced by our new potentials as well. 



The theoretical results together with with experimental data 12Ml6| are shown in Fig. [T] and 
Fig.[2]for one- and two-pole sets of potentials, respectively (we did not take into consideration 



data from 



171 ] with huge error bars). The potentials within each set provide slightly different 



cross-sections, which results in a band instead of a line in the figures. In the same way as 
in j^l we defined "total elastic" K~p cross-section as an integral of the differential cross- 
section over the —1 < cos 6? < 0.966 region due to the singularity of the pure Coulomb 
transition matrix in forward direction. 

Characteristics of kaonic hydrogen atom (AEu, Ti s ) provided by our potentials are situ- 



ated within la KEK 
DEAR collaboration 



18] experimental region, see Fig. [3J Experimental results obtained by 
19 1 are also shown. See J44] as well. 



Typical resonance behaviour manifests itself in 7r°S° elastic cross- sect ions, corresponding 



to one- and two-pole sets of KN — 7rS potential, see Fig. H] (bands, consistin g o: 



individual 



lines). All resonance maxima are situated near Particle data group (PDG) 22| values for 
the A (1405) resonance mass and width 



m pdg 

iW A(1405) 



1406.5 ±4.0 MeV, 



V PDG 
1 A(1405) 



50 ± 2.0 MeV. 



(18) 



PDG mass of the A(1405) resonance is also shown at the figures. Strong pole positions and 
widths are slightly different for the potentials within one- and two-pole sets of potentials. 



TABLE I: Parameters of the representative one- and two-pole KN — -zrS potentials: range f3 a 
(independent on two-body isospin I), strength A" , and additional parameter s of the two-pole 
model. 









(3^ 


\KK 


A 


\ 7T7T 

A 


Ap 




A 7T7T 


s 


Trone- 
KN 


-pole 


3.41 


1.62 


-1.2769 


0.5586 


0.2024 


1.0623 -^0.3251 


1.8315 


1.7158 


0.0000 


x j-twO 

KN- 


-pole 
-ttE 


3.72 


1.00 


-1.6588 


0.4672 


0.0072 


0.7329 -^0.2967 


1.5357 


1.0744 


-0.8433 



For a more detailed description of the properties of our models of KN — 7rS interac- 
tion we chose two "representative" potentials: one with one-pole and another with two-pole 
structure of A(1405) resonance. Parameters of the potentials are shown in Table [U the cor- 
responding observables - in Table [TTJ The latter contains the above mentioned 7, it^s ratios 
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and kaonic hydrogen characteristics AEi s , r± s together with positions of the strong poles 
z\ and z%. Potentials with equal z\ values were chosen as representative ones. The strong 
K~p scattering lengths a K - p , exactly corresponding to the kaonic hydrogen observables, are 
shown in Table [III Since both or— P and (Ai? ls ,r ls ) were obtained by exact solution of 
Lippmann-Schwinger equation, the relation between them does not correspond to any com- 
monly used approximate formula. From Figs. [TJ El [3] and Table [III it is seen, that our new 
one- and two-pole KN — 7rS potentials reproduce all experimental data within experimental 
errors indistinguishably well in the same way as the ones in jsj . Therefore it is not possible 
to give preference to one of the two versions. 

TABLE II: Physical characteristics of the representative one-pole and two-pole potentials (full 
version with physical masses): strong pole(s) position(s) z\ (and Z2), kaonic hydrogen Is level shift 
AEis and width Ti s , K~p scattering length a^--„, and threshold branching ratios 7 and R^t,- 





t 7-one— pole 


■rytwo— pole 
KN-ttY: 


z x (MeV) 


1409 - i36 


1409 - i36 


z 2 (MeV) 




1381 - il05 


AE ls (eV) 


-377 


-373 


r ls (eV) 


434 


514 


a K - p (fm) 


-1.00 + i0.68 


-0.96 + i0.80 


7 


2.36 


2.36 




0.709 


0.709 



In addition, we checked several arguments, which were presented in support to the idea 
of the two-pole structure of A(1405) resonance. One of them is a difference between charged 
7rS cross-sections, which is seen in different experiments, such as CLAS {23)]. In order to 
check the assumption, that the difference is caused by the two-pole structure, we plotted 
7r + £~, 7r~S + , and 7r°S° elastic cross-sections. The result is seen in Fig. the cross-sections 
are different and their maxima are shifted one from another for both one- and two-pole 
versions of KN — 7rS potential. Therefore, it is not a proof of the two-pole structure, but a 
manifestation of an isospin-breaking effect and different background. 
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Another argument for two-pole structure comes from the fact, that the poles in a two- 
pole model are coupled to different channels. Indeed, gradually switching off the coupling 
between KN and 7r£ channels turns the highest of the poles into a real bound state in KN, 
while the lowest one become a resonance in uncoupled 7rS channel. Such a behaviour was 
demonstrated in several papers, see e.g. Accordingly, it was suggested, that the poles 
of a two-body model manifest themselves in different reactions, in particular, KN — KN, 
KN — 7rS, and 7rS — 7r£ amplitudes should "feel" only one of the two poles. We checked 
the hypothesis, the results are demonstrated in Fig. [6j Indeed, real parts of KN — KN, 
KN — 7rS, and 7rS — 7r£ amplitudes in / = state cross real axis at different energies, but 
it is true for both versions of the potential. In fact, the difference is even larger for the one- 
pole amplitudes. In our opinion, the effect is caused by different background contributions 
in the reactions independently of the number of poles. Consequently, a proof of the two-pole 
structure of the KN — 7rS interaction does not exist. 

Coulomb interaction, directly included into two-body Lippmann-Schwinger equations, 
was necessary for reproducing experimental data on kaonic hydrogen atom. However, in the 
K~d scattering length calculations it is expected to play a minor role and can be omitted. 
We also neglected the difference between physical masses in isodoublets for K~d system. 
The physical characteristics of KN — 7rS system, calculated with isospin-averaged masses 
for K and iV using the obtained sets of A" , /3f parameters are shown in Fig. [3l Fig. [3, and 
Table ED 

TABLE III: The same as in Table HT1 but with averaged masses of the particles. 





t rone- pole 


irtwo- pole 
KN--KY. 


zi (MeV) 


1409 - i36 


1409 - i36 


z 2 (MeV) 




1381 - il05 


AE ls (eV) 


-316 


-295 


T u (eV) 


414 


491 


a K - p (fm) 


-0.80 + i0.62 


-0.72 + i0.73 


7 


4.18 


4.54 


R-kT, 


0.761 


0.768 
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"Averaged" points (Ai?i s ,ri s ) for the one- and two-pole representative potentials in 
Fig. [3] are shifted to the smaller |Ai£ ls | values relative to the "physical" ones. However, 
they remain inside la KEK region. Fig. [7] demonstrates "averaged" and "physical" cross- 
sections, where "averaged" ones naturally do not show threshold behaviour at laboratory 
momentum Pi a b, corresponding to K°n threshold. However, differences between "averaged" 
and "physical" cross-sections are visible only in the near-threshold region, where there is no 
reliable experimental data. 

Finally, we see by comparing Table II I II with Table [Til that strong pole positions remain 
almost unchanged. Scattering length a^- p changes for both versions of the potential mainly 
due to the confluence of the K~p and K°n thresholds into one KN threshold. Accord- 
ingly, threshold branching ratios 7 ( TT5"j) and R n ^ ffTB"]) are the only observables, which are 
considerably changed after introducing isospin-averaged masses instead of physical ones. 



B. Nucleon-nucleon potentials 

Antisymmetrized three-body equations for K~d system with s-wave interactions contain 
only spin-triplet NN interaction. We used different NN potentials in order to investigate 
dependence of the K~d scattering length on nucleon-nucleon interaction models. One of 



them is a two-term separable AW potential 24j], which reproduces Argonne V^18 25j phase 



shifts and, therefore, is repulsive at short distances. The potential, which will be called 
TSA, is described by Eq. f[T3"j) with N a = 2 and a — /3 — 1 (the AW interaction is obviously 
is diagonal in particle indices). Two versions of the potential (TSA-A and TSA-B) with 
slightly different form-factors were used: 

»(tr W = E tf > fcr(m) = l,2 (19) 

n=l ^^(m)fi/ 
3 B 2 B 

B,NN / 7 \ _ V" T(l)n n B ' NN (lA - V 7(2) " 

TSA-A and TSA-B potentials in the 3 Si state yield the following scattering lengths and 
effective radii 

a A (np) = -5.402 fm, r^(np) = 1.754 fin, (20) 
a B (np) = -5.413 fm, r* s (np) = 1.760 fm, (21) 
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and give correct binding energy of the deuteron E^cu = —2.2246 MeV. 



We also used one-term PEST potential (Eq. ffl3l) with N a = 1) from Ref. 26], which is 
a separabelization of the Paris model of A A interaction. The strength parameter of PEST 
A = — 1, the form-factor is defined by 

1 6 c NN 

rt) = 2A^H^7' (22) 



where the constants and j3^f are listed in Ref. [26j. PEST is equivalent to the Paris 
potential on and off energy shell up to Ei^ ~ 50 MeV. It reproduces the deuteron binding 
energy E^ eu = —2.2249 MeV, as well as the triplet and singlet AA scattering lengths, 
a( 3 Si) = —5.422 fm and a( 1 S ) = 17.534 fm, respectively. 

The 3 Si phase shifts for the three A A potentials are shown in Fig. [8] together with 
characteristics of the Argonne V18 model. Almost indistinguishable lines correspond to the 
two-term TSA-A and TSA-B potentials, which are very good at reproducing the Argonne 
V18 phase shifts. Their crossing of the real axis is a consequence of repulsion at short 
distances. One-term PEST A A potential does not have such a property, but at low energies 
its phase shifts are also close to the "etalon" ones. Therefore, the two-term A A potentials, 
increasing the number of equations in a three-body system, reproduce properties of NN 
interaction better than the one-term potential. 



C. SiV - AiV interaction 

The £A interaction, which is coupled with AA channel in I = | isospin state, is usually 
assumed to be spin-dependent [27|, [28(. We constructed new versions of EA — AiV potential 
in such a way, that it reproduces existing experimental data 29l-l33|. One-term separable 
potentials, described by Eq. (ll3j) with N a = 1 and Yamaguchi form-factors 

1 



(23) 



were used for the two possible isospin states. But the number of channels is different for 
I — k and / = |. The parameters of the one-channel (a = B = 1) UN interaction with 



isospin I — \ were fitted to the S + p — > S + p cross-sections. Isospin one-half EA is coupled 
to AA channel, therefore, at first we constructed a coupled-channel potential of the I = \ 
SA — AA interaction. The channel indices a, j3 — 1, 2 in ffl3|) denote SA and AA channel, 



12 



correspondingly. The coupled- channel I = \ potential together with the one-channel I = | 
potential reproduces the E~p — » E~p, S~p — y E°n, £~£> — >• An, and Ap —y Ap cross-sections. 

TABLE IV: Range f3 a (independent on two-body isospin I) and strength X"^ parameters of the 
two T,N — AN potentials: V^ dcp and y Sind (g stands for the spin). Scattering lengths af of Y<N 
and A A systems are also shown (in fm). 







2 


pKN 


A? E 

2 


A? A 

2 


A AA 

2 


2 


2 


2 


t rSdep 
V S=0 


1.25 


-0.0244 


0.62 


-1.9956 


1.1408 


-0.7148 


-1.90 + i 0.08 


3.18 


1.26 


t rSdep 
V S=1 


0.50 


-0.0007 


1.03 


-0.0008 


0.0185 


0.0000 


-3.17 + i 1.30 


1.63 


1.57 


^Sind 


0.74 


-0.0032 


0.74 


-0.0011 


0.0254 


0.0190 


-2.40 + i 0.85 


1.95 


-1.47 



Two versions of I = \ EN - AN and I = § EN potentials were constructed: one is spin 
dependent l/ Sde P ; the other V d is independent of spin. Both perfectly reproduce all existing 



experimental data 29H33j on EAf and AAf cross-sections, as is seen in Fig. [9j Parameters of 
the potentials are shown in Table IIVI together with scattering length values and 



a _ 

2 



AA The scattering lengths of spin dependent potential y Sde P are in qualitative agreement 



with those provided by more complicated models of EN interaction 27J, |28j. The only 
exception is af^ with S = 1, having opposite sign, which, however, is the same as that 
given in previous versions of the same advanced potentials (our definition of the sign of 
a scattering length is opposite to those, used in the mentioned articles). The scattering 
lengths, of the spin independent potential 1/ Smd 5 are not in such a good agreement, but 
having in mind, that the scattering length is not a directly measurable quantity, we do not 
consider this difference as a serious defect. 

For the three-body K~d calculations, however, where a channel containing A is not 
included directly, we need not a coupled-channel, but a one-channel EN model of interaction 
in I = | state. Due to this, we additionally constructed an exact optical V SN '° pt and simple 
complex T/Siv,Compiex potentials, corresponding to the obtained I — \ EN —AN potential. As 
was discussed at the beginning of the Section IHIj. the exact optical potential has an energy 
dependent strength parameter defined by Eq. ffT4l and exactly reproduces the elastic EA^ 
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amplitude of the corresponding two-channel potential. Parameters of the y SAr > Com P lex were 
found in such a way, that the simple complex potential gives the same scattering lengths, 
as the two-channel potential. Thus, the exact optical and the simple complex XA(— AiV) 
potentials in / = | state and the one-channel EiV potential in / = | were used during three- 
body calculations. The second channel in brackets (—AN) underlines, that the one-channel 
potentials correspond to the coupled-channel one. 



IV. APPROXIMATE METHODS 



Apart from the full coupled-channel calculation, we performed checks of several approx- 
imate methods, usually used for K~d scattering length problem, as well. It is obvious, that 
a comparison between the full and approximate results is meaningful only if it is performed 
with equal two-body input. 



A. One-channel AGS calculations 



In order to investigate the importance of direct inclusion of ttT,N channel we performed 
one-channel AGS calculations as well. It means, that we solved Eq. (11 01) with a — /3 — 1, 
thus, only KN and NN T-matrices enter the equations. We constructed the exact optical 
and two simple complex one-channel KN(— 7rE) potentials approximating the full coupled- 
channel one- and two-pole models of interaction. As mentioned at the beginning of Sec- 
tion [TTll the exact optical potential V KN, ° pt provides exactly the same elastic KN amplitude 
as the coupled- channel model of interaction. Its energy-dependent strength parameters are 
defined by Eq. (fl4T) with a, (3 = 1,2 stands for KN and 7rE channels, correspondingly. 

For the simple complex potentials we used range parameters j3 KN of the coupled- channel 
KN — 7r£ models of interaction. The complex / \i 1 ' Com P lex constants were obtained in two 
ways. The first version of the simple complex KN potential y^^< Com P lex reproduces K~p 
scattering length ax- P and pole position Z\ of the corresponding coupled-channel version of 
the potential, shown in Table UTTl The second one, ^(f'^)' Complcx provides the same 1 = and 
1 = 1 isospin KN scattering lengths as the full KN — 7rS: 

4 n ^J = -1.60 + i 0.67 fm, a™-J2o = "1-62 + % 0.78 fm, (24) 
4T/=i e = -°- 004 + 1 °- 57 fm ' a ST=°l e = 0-18 + z 0.68 fm. (25) 
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B. Fixed center approximation 



So-called "Fixed center approximation to Faddeev equations" (FCA) introduced in [34 1 
is a variant of FSA or a two-center formula. Fixed-scatterer approximation (FSA) or a two- 
center problem assumes, that the scattering of a projectile particle takes place on two much 
heavier target particles separated by a fixed distance. The motion of the heavy particles is 
subsequently taken into account by averaging the obtained projectile-target amplitude over 
the bound state wave function of the target. Therefore, the approximation is well-known 
and works properly in atomic physics, were an electron is really much lighter than a nucleon 
or an ion. Having in mind, that the antikaon mass is only twice smaller than the mass 
of a nucleon, we can expect, that FSA hardly can be a good approximation for the K~d 
scattering length calculation. 

I — I 

The FCA formula was obtained in [34J from Faddeev equations in a very strange way. 
Proper derivation of a FSA formula starting from the same equations was done much earlier 



m 



35| . it can also be found in [36] together with several versions of the FSA formula. 
Fixed scatterer approximation for the calculation of a K - d scattering length using separable 
potentials was used in [37]. 

A novelty of the FCA formula of 34J is introduction of "isospin breaking terms" , which, 
according to the authors, come from K°n two-body particle channel introduced in addition 
to K~p. However, the inclusion of the K°n channel is questionable, since, as it was already 
mentioned in Section [Til an terms, connected with this subsystem, automatically drop out 
from the Faddeev system of equations after antisymmetrization. 

Necessity to go beyond FCA formula for the K~d system was clearly demonstrated in 



were the unstable character of the FCA results was pointed out. However, the formula is 
still being used, for example in 39[, that is why we decided to check the approximation. 
We used the same two-body input as in AGS equations in order to make the comparison as 
adequate as possible. 

First of all, we used the scattering lengths provided by our coupled- channel KN — 7rS 
potentials and the deuteron wave function corresponding to our TSA-B iViV potential in 
FCA formula Eq.(24) from 34J. Secondly, all K°n parts were removed from the formula 
because they do not enter AGS equations. Finally, we took into account the fact, that the 
FCA formula was obtained for a local KN potential, while separable KN — 7r£ potentials 
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were used in our Faddeev equations. The corresponding changes in the FCA formula were 
made 45[. Therefore, the two-body input for the FCA formula was equivalent to the input 
for the AGS calculation. 



V. RESULTS AND DISCUSSION 

The results of the full coupled-channel calculations of the K~d scattering length using sets 
of one- (empty circles) and two-pole (empty squares) versions of the coupled-channel KN — 
7rE potentials are shown in FigJTUl The calculations were performed with V^ A ~ B and exact 
optical V^ p, ° pt . The K~p scattering length values obtained with the two representative 
KN — 7rE potentials are: 

4 n -7° le = -1-49 + z 0.98 fm, (26) 
a K-7° le = -1-57 + i 1.11 fm. (27) 

Results of previous Faddeev calculations of the same system (filled squares) together with 
two FCA results (crossed squares) are also shown in the figure. 

It is seen, that while two-body data do not allow to distinguish between one- and two-pole 
versions of "the main" KN — 7rE interaction, the three-body sets of results differ sufficiently 
for such a task. Therefore, in principle, it would be possible to favor one version of the 
KN — 7rE potential by comparing with an experimental result. However, direct measure- 
ment of K~d scattering length is impossible. Moreover, it is not absolutely clear, whether 
the difference between the two sets of the ax-d results is much more than theoretical uncer- 
tainties, caused mainly by the uncertainties of the KN model of interaction. In any case, a 
calculation of Is level shift A£?deu,is an d width Tdeu.is of kaonic deuterium atom, correspond- 
ing to the obtained cix-d values, is necessary for comparing with experimentally measurable 
values. The parameters of kaonic deuterium were being measured by SIDDHARTA ex- 
periment, unfortunately, without any results. Due to this our next step will be making 
predictions for Ai?deu,is and rd eu ,is observables. 

As is seen in Fig [TOl our ax-d results are close to the other K~d scattering lengths 



obtained in 



40| and 



41] within coupled-channel Faddeev approach. The result of 42], 
obtained by a one-channel Faddeev calculation with a zero-range one-channel KN potential, 
has much smaller absolute value of the real part than all other a K - d . On the contrary, the 
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authors of 



38j, who performed Faddeev calculations using NN interaction with <i-wave 



component, obtained K d scattering length with quite larger absolute values of both real 
and imaginary parts. 



The ttK-d value of 



34j . significantly different from all others, was calculated using the 



FCA formula, obtained in the same paper and already discussed in Section HVBl We chose 
the result calculated in isospin basis for the comparison. One more paper, where FCA 



formula was used, is [39] . where, however, the result was obtained by simply applying of two 
approximate formulae. The second one is the corrected Deser formula, used for calculation 
of the KN scattering lengths, entering the FCA. One of the representative a^-d values 
from [39(|, having the largest possible imaginary part, is shown in FigfTUl 

It is hard to compare all a K - d results, because the methods of treatment of the three-body 
problem and two-body inputs are different in the mentioned works. In order to investigate 
separate effects of several approximations we performed approximate calculations, as de- 
scribed in Section [iVl The obtained one-channel AGS ax-d values (Section IIV Al) together 
with FCA (Section IIV Bl) and the representative coupled-channel AGS results of the K~d 
scattering length calculations are shown in Fig. []J]for one- and two-pole versions of KN — tiH 
interaction. It is important, that all results in the figure were obtained with equivalent two- 
body input, including the neglect of isospin-breaking parts in the original FCA formula. 

It is seen from FigHTl that all approximations are more accurate for the one-pole version 
of the KN interaction than for the two-pole variant. The one-channel AGS calculation with 
exact optical KN potential (empty symbols), giving exactly the same KN — KN amplitude 
as the corresponding coupled-channel potential, turns out to be the best approximation. The 
result, obtained with the simple complex potential V^^ ,complex (vertically crossed symbols), 
reproducing 1 = and 1 = 1 KN scattering lengths, underestimates the absolute value of 
the real part of the a^-di especially for the two-pole version of KN interaction, but have 
rather accurate imaginary part of the K~d scattering length. Another one-channel AGS 
calculation with simple complex potential y^^> com P lex ; reproducing K~p scattering length 
and pole position zi, gives rather inaccurate result (half-empty symbols) as compared with 
the coupled-channel AGS values (filled symbols). 

The scattering length ax-d of [42] was obtained from one-channel Faddeev equations 
with a complex potential. However, the underestimation of the absolute values of its real 
part in comparison to other Faddeev calculations is so large, that it cannot be explained 
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by the method only. The most likely reason of the difference lyes in the properties of the 

KN potential used in |42| . First of all, the potential has a very large mass of the K~p 

quasibound state (1439 MeV), which, therefore, is situated above the K~p threshold. In 

addition, the width of the state (127 MeV) is much larger than the PDG value (50 MeV) as 

well as the width of our z\ pole (72 MeV). 

It is hard to understand the results obtained in [33]. While all formulae are written 

for many-channel Faddeev equations, the most of the dependences and even "the best" 

ciK-d value were obtained within a one-channel Faddeev calculation including KN and NN 

interactions only. Since the elastic part of the coupled-channel KN T-matrix was used, the 

result is equivalent to a one-channel Faddeev calculation with an exact optical potential. 

I 

But even the full coupled-channel calculation was performed in [38] with non-unitary KN 
T-matrices OSA and OS1, since channels with 77-mesons, entering the two-body T, were 
omitted in the three-body equations. It is not clear, why the one-channel calculations of 38 ] 
give so large difference in imaginary parts of ax-d obtained in isospin and particle basis and 
with and without (i-wave in NN. The result of a calculation, in principle, should not depenc 
on the chosen basis, in addition, the very recent results of K~d scattering calculations {43] 
demonstrated, that inclusion of physical masses into Faddeev equations change ax-d by 
several percents only. 

The results of using the FCA formula without isospin-breaking effects (crossed symbols) 
stays far away from the full calculation, as is seen in FigJTTl While errors for the imaginary 
part are not so large, the module of the real part is underestimated by about 30%. There- 
fore, our calculations show, that FCA is a poor approximation for K~d scattering length 
calculation, and the accuracy is lower for the two-pole KN model of interaction (the most of 
chirally-based models of KN interaction have two-pole structure). Even the original FCA 
formula does not give such a large K~d scattering length as a^-d from 34j, which, therefore, 
is caused by too large input KN scattering lengths. As for the values of 39], they suffer 
from cumulative errors of two approximations and using of DEAR results on kaonic hydro- 
gen characteristics. As was shown in 3|, the error of the corrected Deser formula makes 
about 10%, while the problems with DEAR experimental data were also demonstrated in 
the paper and in other theoretical works. 

We investigated dependence of the full coupled-channel results on NN and £iV(— AiV) 
interactions as well. The dependence of a^-d 011 nucleon-nucleon interaction is demonstrated 
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in Fig [121 were the results obtained with TSA-A, TSA-B, and PEST V , are shown. The 
representative sets of one- and two-pole KN — 7rS potentials were used together with exact 
optical V Sdep ' 0pt EN (-AN) model of interaction. We see from the Figure, that the difference 
is very small even for the potentials with and without repulsion at short distances (TSA and 
PEST, correspondingly). Therefore, the s-wave iViV interaction, which is used in the present 
calculation, plays minor role in the calculation. Most likely, it is caused by relative weakness 
of the NN interaction as compared to KN from the viewpoint of a much deeper quasibound 
state in the latter system (Ej <N m —23 MeV for our potentials) than the deuteron bound 
state (-Edcu ~ — 2 MeV). We do not expect much larger effect from higher partial waves in 
NN as well. 

We also looked at the dependence of a K ~ d on T,N(—AN) interaction. The K~d scat- 
tering lengths obtained with the exact optical and the simple complex versions of the spin 
dependent V Sdep and spin independent V d potentials are shown in Fig. [131 The repre- 
sentative sets of one- and two-pole KN — nil potentials were used together with TSA-B 
NN potential. The results of the two versions of £iV(— AiV) potential V Sdep and V d in 
exact optical form are very close, while their simple complex versions are slightly different. 
However, the largest error does not exceed 3%, therefore, the dependence of the a>K- p on 
EiV — (AiV) interaction is also weak. 

VI. CONCLUSIONS 

To conclude, we performed calculations of the K~d scattering length using newly ob- 
tained coupled-channel KN — 7r£ potentials with one- and two-pole versions of the A(1405) 
resonance. Faddeev-type AGS equations were used for description of KNN — 7r£iV system. 
We also constructed new coupled- channel UN— AN potentials together with its exact optical 
and simple complex SV(— AiV) versions. Different models of the A^A" interaction: TSA-A, 
TSA-B, and PEST were used. All two-body interactions are described by s-wave separable 
potentials. We investigated dependence of ax-d on NN and SV(— AiV) interaction and 
found, that both dependences are weak. 

We found, that the two sets of the results, obtained with one- and two-pole models of 
A(1405) resonance, are clearly separated one from another, in principle, allowing to give 
preference to one of the KN — 7rS interaction models. However, the question, whether 
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theoretical uncertainties are not of the same order as the differences between the two obtained 
sets of ax-di remains open. In any case, it is necessary to calculate level shifts and widths 
of kaonic deuterium atom, corresponding to the obtained K~d scattering lengths, which can 
be measured, say, by SIDDHARTA-2 experiment. It is assumed to be done in a next paper. 

Among approximate results, the one-channel AGS calculation with exact optical 
KN{—irY) potential gives the best approximation to the full coupled-channel result. On the 
contrary, FCA was shown to be the least accurate approximation, especially in reproduction 
of the real part of the K~d scattering length. All approximations are less accurate for the 
two-pole model of KN — 7rS interaction. 
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FIG. 1: Comparison of the elastic and inelastic K p cross-sections (filled circles) for the one-pole 



sets of the KN — 7rS potential with experimental data 
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16] (data points). The theoretical bands 



are formed by all lines obtained with individual potentials within the set. 
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FIG. 2: Comparison of the elastic and inelastic K p cross-sections (filled circles) for the two-pole 



sets of the KN — 7rS potential with experimental data 
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16] (data points). The theoretical bands 



are formed by all lines obtained with individual potentials within the set. 
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FIG. 7: Comparison of the theoretical K p cross-sections for the representative one- and two-pole 
KN — ttTi potentials with experimental data 
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161 ] (data points). The results obtained with the 



physical masses (solid, dash-dotted line) and the averaged masses (dashed and dotted line) are 
presented. 
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FIG. 11: Comparison of the full results of the a^-d calculations with representative one- (circles) 
and two-pole (squares) versions of KN — ttE potential with approximate results. Values obtained 
with coupled-channel AGS equations (filled symbol), one-channel AGS with exact optical KN 
potential V KN , opt (empty symbol), one-channel AGS with complex V , complex^) (half-empty 
symbol) and V , complex^) (vertically crossed symbol) results are shown. Results of corrected 
FCA formula using (crossed symbol) are also shown. 
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FIG. 12: Dependence of the full if _ c2 scattering length results on NN interaction: a^-d obtained 
with with TSA-A (empty symbol), TSA-B (filled symbol), and PEST (vertically crossed sym- 
bol) NN potentials. One- (circles) and two-pole (squares) representative versions of KN — 7rS 
interaction were used. 
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FIG. 13: Dependence of the full K~d scattering length results on SiV— (AN) interaction. The a^-d 
values obtained with exact optical yS^P-op* (filled square), simple complex T/Sdep,com P iex ( emp ty 
square), exact optical V Smd,opt (half-empty square), and simple complex T/Smd,compicx ( ver tically 
crossed square) are shown. One- (circles) and two-pole (squares) representative versions of KN—irY, 
interaction were used. 
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